Astronomy & Astrophysics manuscript no. cooling 


© ESO 2011 


December 19, 2011 





Catastrophic cooling and cessation of heating in the solar corona 



H. Peter, S. Bingert, S. Kamio 



o 

(N 
o 

CD 

Q 

in 



Oh 
O 



> 

m 



X 



Max-Planck-Institut fur Sonnensystemforschung, 37191 Katlenburg-Lindau, Germany, email: peter@mps.mpg.de 
Received 16 August 2011 / Accepted 12 November 2011 

ABSTRACT 

Context. Condensations in the more than 10'' K hot corona of the Sun are commonly observed in the extreme ultravi- 
olet (EUV). While their contribution to the total solar EUV radiation is still a matter of debate, these condensations 
certainly provide a valuable tool for studying the dynamic response of the corona to the heating processes. 
Aims. We investigate different distributions of energy input in time and space to investigate which process is most 
relevant for understanding these coronal condensations. 

Methods. For a comparison to observations we synthesize EUV emission from a time-dependent, one-dimensional model 
for coronal loops, where we employ two heating scenarios: simply shutting down the heating and a model where the 
heating is very concentrated at the loop footpoints, while keeping the total heat input constant. 

Results. The heating off/on model does not lead to significant EUV count rates that one observes with SDO/AIA. In 
contrast, the concentration of the heating near the footpoints leads to thermal non-equilibrium near the loop top result- 
ing in the well-known catastrophic cooling. This process gives a good match to observations of coronal condensations. 
Conclusions. This shows that the corona needs a steady supply of energy to support the coronal plasma, even during 
coronal condensations. Otherwise the corona would drain very fast, too fast to even form a condensation. 

Key words. Sun: corona — Sun: UV radiation — Sun: X-rays, gamma rays — Sun: activity — Hydrodynamics 



1. Introduction 

Loops in the corona still present us with numerous open 
questions. Besides the open fundamental problem of the 
coronal heating mechanism, it is not clear yet which pro- 
cesses can provide the high densities observed high in the 
corona, well above what would be expected by hyd r ostati c 
equilibrium (e.g. iPeresl Il997t lAschwanden et all . l200ll) . 
One candidate for this is the thermal non-equilibrium that 
leads to catastrophic cooling. This process is initiated if the 
heating is very concentrated at the loop footpoints, and 
the heating at the loop top can no longer balance the losses 
from heat conduction and radiation: the loops starts cool- 
ing, radiatio n becomes more efficien t, and a runaway pro- 
cess sets in (|Kuin fc Martensl . |1982[ ). This process is con- 
sidered to be important in th e formation of prominences 
(|Antiochos fc Klimchukl . Il99ll: iKaroen et all [2006) . The 
role of this catastrophic cooling for normal coronal loops 
has been studied in several numerical experiments that 
employ one-dimensional models for (semi-circular) loops. 
iMuiller et all ()2003l |2004[ ) have investigated the cyclic and 
chaotic behaviour for loops of different lengths and derived 
observable quantities such as extreme ultraviolet emission 
line spectra. They relat e these events to "coronal rain" 
(fMti ller ct al., 2005). An toUn et aD (|2010[ ) proposes that 
"coronal rain" could be a marker of the coronal heating 
mechanism. 



In a recent study iKlimchuk et al.l (|2010D have been 
able to show that possibly the catastrophic cooling asso- 
ciated with the thermal non-equilibrium does not play a 
major role in active region coronal loops. Consequently, 
iKUmchuk erall EoWti argue that the heat input for active 
region loops is probably not concentrated near the loop's 



footpoints alone. The authors furthermore suggest that the 
heating is probably not steady. 

The conclusions of IKlimchuk et al.l ()2010l) on the rel- 
evance of the process of thermal non-equilibrium are still 
under debate, as has become clear in discussions at the re- 
cent Loops- 5 workhop (June 2011). When the loop cross- 
section is allowed to vary, and especially when the heat- 
ing profile is made asymmetric, the nature of catastrophic 
cooling may be signif icantly different than concluded by 
I Klimchuk et al.l (|2010h . It may be premature to rule out 
thermal non-equilibrium for explaining some of the ob- 
served properties of coronal loops (Z. Mikic, priv. comm.). 

Three-dimensional magneto-hydrodynamic models (3D 
MHD) of the corona show a horizontally averaged expo- 
nential drop in the heating rate wit h altitude, with a scale 
height of th e order of 5 Mm (e.g. iGudiksen fc Nordlundl . 
I2OO2I l2005al lb[) , and the heating is fou nd to be intermittent 
in time and space (iBingert fc Pet^ |2011D . The assump- 
tions of the IKlimchuk et al.l i2010) study concerning the 
spatial and temporal distribution of the heat input are con- 
sistent with the average drop of the heating rate found in 
the 3D MHD models. This type of 3D MHD model can be 
considered realistic, since they successfully reproduce nu- 
merous properties observed in imaging and spectroscopy, 
such as the emission measure distribu tion or the aver- 
age transition region Doppl er shifts Ce.g. lPeter et al.ir2004 
I2OO6I : [Zacharias et al.ll20lil ). Because these 3D models pre- 
dict a heating scale height greater than what is needed for 
catastrophic cooling, this is further evidence that extreme 
concentration of the heating rate towards the footpoints is 
not a general feature of coronal heating. 

One can consider catastrophic cooling as a process that 
will happen only in a sub- volume of the corona, where the 
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heating (for some time) is concentrated much more towards 
the loop footpoints than on average. These catastrophic 
cooling events are observed and seem to be a common fea- 
ture. Examples of detailed obs ervationa. 1 studi es of such 
events are the investigations by ISchriiveJ ()2001[ ) in a com- 
parison of filtergr ams in Lyman-g and in the coronal band 
around 171 A, bv Ide Groof et all (|2005l) when comparing 
the 17 1 A band with H-a data, or bv lAntolin fc Verwichtd 
(|2011h on the role of these condensations for loop oscilla- 
tions. 

Most recently, cooling events have been reported by 
iKamio et al](|20lH ) using data from the Atmospheric Image 
Assembly on t he So lar Dynamics Obs ervatory (AIA/SDO , 
iLemen et al.l . |2011D . In particular, iKamio et al.l ()2011[ ) 
present cases with peculiar light curves: the 131 A and 
171 A representing plasma at logT[K]«5.7 and 5.9 show 
narrow single peaks in time. In contrast, the 193 A and 
211 A channels with a maximum contribution for plasma 
at logr[K]R::6.2 and 6.3 show enhanced emission for much 
longer times with a more complex light curve, including 
peaks before and after the peaks in the "cooler" chan- 
nels. A more detailed presentation of this lightcurve and 
the comparison to our model can be found in Sect.[5l The 
lightcurve found by Kamio et al. (2011) is counter-intuitive 
for a situation where the plasma cools, and one would ex- 
pect a brightening in the channels according to their tem- 
perature of maximum contributio n. A simplistic model was 
presented bv lKamio et al.l (|201l[ ) to give a first attempt to 
explain this peculiar behaviour assuming a cooling at con- 
stant pressure. 

In the present study we investigate the cooling of the 
plasma i n order to underst and the peculiar light curves 
found bv lKami~ et al.l ()20lH) . We perform one-dimensional 
loop models starting from an equilibrium model that is hot 
enough that it does not show significant count rates in any 
of the above channels. We then follow two scenarios: (1) a 
catastrophic cooling case where the heating is concentrated 
at the foot points, and (2) a heating off /on case where the 
heating is shut off completely for a while and the turned on 
again. In both cases the plasma will cool, but only the first 
case will be compatible with the observed cooling events. 

After a short outline of the model in Sect. [5] we present 
the hydrodynamic results (Sect.|3]) and the observable AIA 
signatures (Sect.|3]) for both c ases, before we discu ss a com- 
parison to the observations of IKamio etall (poTTh m more 
detail in Sect. [5] and conclude the paper. 

2. ID loop model and synthesizing observations 

To model the dynamics of a coronal loop we solve the 
mass, momentum, and energy balance in a ID model 
for a semi-circular loop with constant cross-section. For 
the numerical experiments w e employ the Pencil code 
(jBrandenburg fc Dobled . l2002f !Fl with modifications to ac- 
count for the physics of the corona (Bingert & Peter, 2011). 
The equations (in their full 3D form) and a brief descrip tion 
of the method can be found in iBingert fc Peteil (120111 ). 

For the loop model we assume that all variables are a 
function of the coordinate s along a magnetic field line, 
which is assumed to be semi-circular. The velocity is only 
along the loop. At both ends we impose boundary condi- 
tions reflecting the photosphere; i.e., we prescribe the den- 

^ pencil-code.googlecode.com/ 



sity and the temperature . The velocity is s e t to z ero at both 
ends. In the 3D-model of' Bingert fc Peteil (|201l[ ). the heat- 
ing through Ohmic dissipation is self-consistently described 
and results in a roughly exponential decrease in the heating 
rate in the coronal part with a scale height of about 5 Mm. 
In this ID model we simply define the heating rate to be ex- 
ponentially decreasing and constant in time (see Sect. 12.1]) . 
The energy equat ion also accounts f or optically thin radia- 
tive losses (fro m ICook et all . |1989[ ) and heat conduction 
(jSpitzed . Il962[ ). which is essential to properly model the 
thermal non-equilibrium leading to catastrophic cooling. In 
the momentum equation, gravity is taken into account. 

The model produces a hot corona above a cool pho- 
tosphere and chromosphere at both footpoints with a thin 
transition region. The lower dense and cool part of the loop 
is only used as a reservoir for mass and energy, so we do not 
solve the radiative transfer problem, among others. For the 
simulations shown in this manuscript we used 2048 grid- 
points along the loop. 

2.1. Heating rate and equilibrium model 

For the heating rate we prescribe an exponentially decreas- 
ing volumetric heating rate, 

g = goexp(-|) , (1) 

with the heating rate at the lower boundary Qo and the 
scale length of the heating rate A. This is a function of 
the height z and not of the arc length along the loop. The 
energy flux density into the corona -Fh is the flux at the 
base of the corona, 

Fh = / Qdz , (2) 

where ztop is the height of the loop apex. We define the 
base of corona, ^basc, as the height where the transition re- 
gion starts, i.e., where a temperature of 10^ K is reached. 
Essentially, there are two free parameters for the heating 
rate, namely the energy fiux into the corona, .Fh, and the 
degree of concentration of the heating towards the foot- 
points parameterized through the heating scale length A. 

To obtain a thermal non-equilibrium causing a catas- 
trophic cooling event, the heating scale length A has to 
be short enough compared to the loop length. Because we 
want to compare an event with another one without catas- 
trophic cooling, we choose A=2 Mm. This provides a sta- 
ble solution (while A=l Mm shows a loss of equilibrium) 
for the loops with a length of L=120Mm that we study 
here. This loop len gth is motivated by the observation of 
IKamio et al.l ()2011|) . who find the condensation at an apex 
height of about 40 Mm. 

Because we intend to investigate a loop that is initially 
too hot to be seen in the SDO/AIA 193 A or 211 A channels, 

the apex temperature should be at least logr[K] w 6.5. 
This automatically sets the required energy fiux to heat 
the corona. In our equilibrium (numerical) model we choose 
a heating rate resulting in a coronal base at Zbaso ~ 
3.4 Mm with a rather high coronal energy input of -Fh ~ 
15 000Wm~^. This gives the required apex temperature 
of logT [K] w 6.5. This choice of parameters automatically 
sets the pressure in the loop and thus the density at the 
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apex, which is about logri [cm ■^J « 9.5 in the equilibrium ^ 
model. .'>< 
The heat input (and the loop length) fully determines 7 
the temperature and the den sity. This is also refl ected by 
the scaling laws as derived bv lRosner et al.l (|1978[ ). 



f [K] = 1700 (FH[Wm-2])2/^ {L[m]f^^ 



[cm-3] = 3.9-101° (FH[Wm-2])4/^ (i [m 



-3/7 



(3) 
(4) 



In this modi f ied fo rm the (constant) heating rate E^i of 
iRosner et al.l (|1978[ ) is translated into the energy flux den- 
sity into the corona, Fh = En {L/2). For the above values of 
energy flux and loop length used for our equilibrium model, 
these scaling laws give comparable values for the tempera- 
ture and density at the loop apex as the numerical model, 
as should be expected. 



2.2. Numerical experiments 

We started the two different numerical experiments from 
an equilibrium model with the parameters for the heating 
rate as mentioned in Sect. 12. l1 In particular, the heating 
scale length is A = 2 Mm. We let this model run for a long 
time to be sure that it represents an equilibrium solution. 
"Long" means much longer than the coronal cooling time, 
which is around hour. 

The first systematic study for loops heated with the 
exponenti ally decreas i ng he ating rate used here was per- 
formed bv lSerio et al.l (jl98lfl . For small loops (shorter than 
the pressure scale height), they could find solutions for 
their static model with a temperature peak at the apex 
only if the heating scale height was more than half the 
loop length. For longer loops they found this threshold is 
a third of the pressure scale height. Thus this would be 
about 50 Mm in o ur case (of long loops), and according to 
ISerio et al.l (p8ll) we should not find a stable solution for 
A = 2 Mm. However, later time- depend ent models such as 
the ones by iMiiller et all (|2003l I2004D have found stable 
solutions, also at lower values for A. Whether the thermal 
non-equilibrium sets in depends on the balance of heating 
and radiative losses (at the loop apex). This is a complex 
non-linear balance, because decreasing the scale length, and 
thus decreasing the heating at the apex, also leads to a 
lower density at the apex, hence less radiation. Therefore 
the quantitative results for the thre s hold value of A based 
on the static models of ISerio et al.l ()1981[ ) do not hold in 
general. Furthermore, for the present study it is not so 
important on which heating scale height the thermal non- 
equilibrium sets in, but only that it sets in at all. 



2.2.1. Catastrophic cooling 

In one experiment we decrease the heating scale length to 
A = 1, while keeping the energy flux into the corona, Fh 
as deflned in ([2]), constant. Thus we changed Qo in ^ 
accordingly. All other parameters are kept exactly the same. 

This setup will result in a thermal non-equilibrium near 
the top of the loop, leading to catastrophic cooling and the 
formation of a condensation. After a while the loop recovers 
its former temperature (and density). This happens in a 
similar fashion to the one described bv IMiiller etal] ()2003l 
12004 . 



2; 
Q 




5.5 6.0 6.5 
temperature log T [K] 

Fig. 1. Temperature response functions of the AIA chan- 
nels as stored in Sol arSoft at the t ime o f publication of 
this study; following iBoerner et al.l (|201lD . Noted are the 
centre wavelengths in A and the peak of the temperature 
response in logT[K]. The secondary peaks at low temper- 
atures in the 193 A and 211 A channels are indicated by 
arrows. See Sect. 12. 31 

2.2.2. Heating ofF/on 

In the other experiment, we kept the same heating scale 
height, but practically shut off the heating by decreasing 
the heat input by many orders of magnitude (10^^) from one 
time step in the simulation to the next. All other param- 
eters were kept as in the equilibrium model. Starting with 
the equilibrium solution this implies that the loop cools 
down and the plasma will drain. 

After some time we resumed the original heating rate, 
and the loop found its former equilibrium. We waited before 
turning the heating on for a time comparable to the time 
it takes the catastrophically cooling loop fSect. ETTTl) to 
recover from its cooling event. 



2.3. AIA synthetic data 

To allow for a comparison with SDO/AIA observations, we 
synthesized the count rates to be expected if AIA were to 
observe the structure we model. For this we use the temper- 
ature response functions as provided in the SolarSoft pack- 
agfl The temperature response function or kernel Ki{T) 
is basically the integral of the contribution function of the 
lines in the respective bandpass i weighted with the effective 
area as a function of wavelength. The ke rnels Ki{T) are de- 
scribed i n detail in Boerner et al.l (|201ll ) and are displayed 
in Fig.[T] iBoerner et al.l (|20 111) use d ifferent abundance s and 
another ionization equilibrium than lCook et al.l ()1989D used 
to derive their radiative loss function, which we employ in 
our model. This inconsistency should not significantly affect 
our results, but should be checked in the future. 

The emissivity (i.e., radiated power per volume) at each 
gridpoint of the loop model is given by K{T), with n and 
T the number density and temperature at that gridpoint. 
To derive the AIA counts in the common digital number 
DNpixel^i s~i, one then has to integrate along the line of 
sight. Because the loop model is one-dimensional, we have 



http:/ /www. lmsal.com/solarsoft/ 
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to assign a diameter to the loop to derive the synthetic 
count rate. For the investigations presented here we assume 
a diameter of 50 km. This choice is motivated by the di- 
ameter of the individual loop strands in the m ulti-stranded 
loop model of iPatsourakos fc KlimchuS ()2006[ ). but the re- 
sults presented here do not depend much on this special 
choice. 



3. Results for ID loop models 

3.1. Loop undergoing catastrophic cooling 

Starting from the equilibrium model, the loop gradually 
evolves after having decreased the scale length A of the 
heating rate (cf. Sect. l2.2."T|) . For approximately the first 1.5 
hours the changes in the plasma parameters, i.e., tempera- 
ture, pressure, and velocity are only very small. Therefore, 
in Fig. [2] (bottom row) we show the evolution of these quan- 
tities for this model only after tps90 min. In this early phase 
the loop reaches a high temperature with a flat peak at 
about logT [K] = 6.5 (see Fig.[3l top left). Initially the de- 
crease in the heating rate puts the loop only slightly out of 
equilibrium. Because the initial peak temperature is high, 
the (radiative) cooling time is long, and thus the slow early 
evolution is to be expect e d. This is consistent with ear- 
lier work bv lMiiller et al.l (pOOS. 2004), who find that the 
time between the condensations range from one to several 
hours, depending on loop length and heating scale length. 
Depending on its density and its density contrast to the 
surrounding corona, the initial loop might be visible in ob- 
servations either as an individual X-ray coronal loop (stable 
for well above an hour) or as part of the diffuse background 
corona. 

Finally, at about t=110min the loss of equilibrium be- 
comes increasingly faster and a cool dense condensation 
forms at about t= 120 min. While cooling down, the pres- 
sure across the condensation is roughly constant, which is 
due to the inflow of plasma from both sides along the loop 
leading to an increase in density. This behaviour can be 
found in a quantitative fashion through the profiles shown 
in the left-hand panel of Fig. [3] It is also clearly evident in 
the space-time plot in the bottom row of Fig.[5J where one 
sees the condensation in the temperature plot, but only 
a comparably weak trace in the pressure plot. Of course, 
there is some change in the pressure during the conden- 
sation, about 0.2 in logj^p, equivalent to a factor of 1.5. 
However, this is small compared to the change in temper- 
ature (or density), which is almost 2.5 in log^g, equivalent 
to a factor of 300. Therefore, the catastrophic cooling can 
be considered to be more or less isobaric. 

After hovering for some time close to the apex, the con- 
densation slides down the loop and hits the photosphere. 
In this model run the condensation happens to be located 
just a bit left of the loop top (when looking at Fig. [2), and 
thus the condensation slid to that side. This is caused by 
slight asymmetries, e.g., by waves travelling back and forth 
through the loop following the change in the heating scale 
length at t—0. All of this behaviour of the condensation 
is consiste nt with previous work on catastrophic cooling in 
loops (e.s. iMiiller elan l200,l [2nnl . 



3.2. Loop with heating shut ofF and on again 

For the case of shutting off the heating, the temperature 
and the pressure of the loop respond immediately, which is 
clearly visible in Fig.[2^,c. The temporal evolution at the 
top of the loop shown in Fig.S] (top right) gives a quantita- 
tive measure, showing that both temperature and the den- 
sity fall off roughly exponentially on a time scale of about 
10 min, which is comparable to the sound crossing time. 
Because the energy losses through radiation and heat con- 
duction are no longer replenished, the plasma cools and con- 
sequently loses its support: the loop quickly drains. This is 
visible in Fig.[2i, which shows the velocity along the loop. 

This draining of mass is a pointed difference to the case 
of catastrophic cooling discussed in the preceding section. 
Once the condensation has formed there, it goes on to ac- 
quire mass, which is visible as upflows in both legs of the 
loop (Fig.[5]i starting at tRil20min). In contrast, downflows 
in both legs can be observed in the case with the heating 
shut off (Fig.[5Jl) . This draining is efficient enough that af- 
ter some 50 min the density dropped by about a factor of 
100, and because of the drop in temperature the pressure 
at the apex dropped by about 10^. This underlines that the 
case of catastrophic cooling can be considered as a basically 
isobaric process, at least its pressure change is smaller than 
in the case of shutting off the heating. This is because when 
basically following ([3]) and (|H) the total heating rate sets the 

pressure, poc-F^^^, and the total heat input Fn remains the 
same when there is of catastrophic cooling (cf. Sect. l2.2.T|) . 
This has an important consequence when discussing the 
signatures that are observable with AIA (Sect. 2]). 

After the loop has cooled and drained we switch the 
heating on again (at t«57min) to return to the value of 
the equilibrium model. This leads to a very fast increase 
in the temperature. Within one minute the temperature at 
the apex rises to the original value of about logT[K]«6.5 
(cf. Fig.|3]). The increased heating rate leads to evaporation 
at the loop footpoints, and the upflows fill the loop to its 
former density within some 15 min (cf. Fig. 21 left panels). 
This very strong change in density is also pointedly differ- 
ent from the catastrophic cooling, where the density in the 
loops outside the condensation region remains more or less 
constant (cf. Fig.|31). 

4. Observational signatures of the cooling events 

Because the AIA bandpasses have contributions from mul- 
tiple temperatures, the expected count rates in a cooling en- 
vironment do not necessarily show a peak when the plasma 
reaches the temperature of maximum contribution. This 
has already been pointed out in our observational study by 
employing a simplistic model for a co oling plasma at con - 
stant pressure (Sect. 3.3 and Fig. 6 of iKamio et al.l . [20111 ). 
The numerical loop model presented here supports and sig- 
nificantly refines our previous arguments. 

In the following we analyse the AIA emission, as can 
be expected at the loop apex when looking horizontally at 
the loop edge-on (i.e., the observer is in the same plane 
as the loop). To capture the condensation at the loop top, 
we concentrate on the part that would be covered by one 
AIA pixel, i.e. «0.6" = 450 km. Because we look at the loop 
apex edge-on, this corresponds to an integration over an arc 
length of about 6 Mm to both sides of the apex. We discuss 
only the 304 A, 131 A, 171 A, 193 A, and 211 A channels 
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Fig. 2. Plasma parameters in numerical loop models as a function of arc length along loop and time. The loop is semi- 
circular with the apex at 60 Mm and the footpoints at and 120 Mm. Shown are temperature, density, pressure, and 
velocity (from left to right) for the two numerical experiments. The top row displays the results where the heating was 
shut off (at t=0) and then turned on again at tss57min (see Sect. 13. 2p . In the bottom row, the results for the case of 
catastrophic cooling are shown (see Sect. 13. ip . where t=0 refers to the time when the heating scale length was reduced. 
The condensation sets in at around tRil20min. In both cases the time axis spans about 70min. Positive velocities (shown 
in blue) are in the direction of increasing arc length. 



of AIA. The 335 A channel has a very broad contribution, 
which is not really useful for studying a condensation, and 
for the 94 A channel it is most likely that the temperature 
response function bas ed on CHIANTI calculations has to 
be severely corrected (jAschwanden fc Boerneil 1201 it) . 



count rate comes about because the pressure is roughly 
constant, and thus the density increases. This is due to the 
weighting effect with the density, because the count rate is 
given through the temperature response weighted by the 
density squared. 



4.1. Catastrophic cooling 

All the AIA channels under investigation show very low 
count rates at the initial stage of the catastrophic cool- 
ing experiment (which is by design, see Sect. 12. ip . In Fig. [3] 
(right panel) we show the temporal evolution of the synthe- 
sized AIA counts starting at about t=90min, because the 
early evolution is very slow (cf. Sect. 13. ip . 

"Hot" channels, 193 A and 211 A. In these we first see a slow 
increase in the count rate. This is because the temperature 
is falling and the broad peak very roughly coincides with 
the time when the apex reached the temperature, where the 
temperature response of the respective channel is at maxi- 
mum (at logr[K]«6.2 and 6.3; indicated by vertical dotted 
lines) . The considerable time lag between reaching the peak 
of the temperature response and reaching the peak of the 



After this broad first peak, a very narrow one can be 
seen m the 193 A and 211 A channels (at i«120 min) when 
the condensation sets in. It occurs when the apex tem- 
perature reaches the value corresponding to the secondary 
major contribution to these channels at logT[K]sa5.4 (cf. 
Fig .m. This second peak is due to ions such as Ovi 
(jO'Dwver et al.l[2010[ ). It is much sharper in time, because 
the temperature falls rapidly when forming the condensa- 
tion. Even though the secondary maximum of the temper- 
ature response is a factor of about 10 lower than the main 
maximum (Fig.[T]) and the source volume is much smaller 
(because the transition region is very narrow), the count 
rate during this second peak is comparable to the first peak 
(in the case shown in Fig. [3] only a factor of 2 lower) . This 
is because the condensation evolves basically at constant 
pressure (cf. Sect. 13. 1| ): when reaching logr[K]w5.4 the tem- 
perature drops by a factor of 10, thus the density increased 
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arc length [Mm] time [ minutes ] 



Fig. 3. Plasma properties and synthesized emission of AIA channels for ID loop model undergoing catastrophic cooling. 
The left column shows snapshots of the profiles of temperature, logT, density, logn, pressure logn, and velocity, v along 
the loop (w>0 in direction of positive arc length). The right column shows the temporal evolution of the temperature 
and density at the loop apex (top right). The other panels in the right column show the expected count rates for several 
AIA/SDO channels if the loop top is observed edge on (one AIA pixel covering apex at 38 Mm height). The short vertical 
markers indicate the time for which the snapshots are shown in the left column. The vertical dotted lines indicate the 
times when the temperature at the loop top (top panel) reaches the temperature of maximum response for the AIA 
channels. These temperatures are given (in logT[K]) with each panel. See Sects. lOI and l4Tl A movie of the temporal 
evolution shown in the left panels is available in the on-line edition. 



by a factor of 10, and the temperature response has to be 
weighted with the density squared! 

Finally, once the condensation has slid away from the 
apex, the plasma resumes its original temperature and den- 
sity, and we see a third maximum in the 193 A and 211 A 
channels, which is basically a mirror of the first maximum. 
Because the density does not immediately fully return to 
its original value, the count rates in this second maximum 
are a bit lower than during the first maximum. 

"Cool" channels, 131 A and 171 A. These cooler channels, 
with a maximum contribution at logT[K]Ri5.7 and 5.9, show 
basically the same as the "hot" channels for the first and the 
third maximums in the count rate. However, the temporal 



evolution in the "cool" channels is faster; i.e., the peaks are 
narrower, because the condensation crosses the maximum 
contribution temperatures faster (with decreasing temper- 
ature the first and the last peaks become increasingly nar- 
row, cf. Fig. [3]). 



The major difference between the "hot" and the "cool" 
channels is that the 131 A and 171 A channels do not show 
the second maximum. This is because these channels show 
a lesser degree of contamination with cooler transition re- 
gion lines; in particular, they do not show any secondary 
peak at temperatures below the main peak (even though 
the contribution shows some extended wing; cf. Fig.[T]). 
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Fig. 4. Same as Fig. [3] but for a loop with heating switched off at t=Omin and on again at sa57min. See Sects. and 
14.21 A movie of the temporal evolution shown in the left panels is available in the on-line edition. 



Hen channel, 304 A. The Heii channel shows a strongly 
asymmetric temporal variation. It shows a first peak when 
the temperature falls below 10^ K (at tRil20min). The sec- 
ond broad peak is from the accumulation of cool material 
near the loop apex as long as the condensation is close to 
the apex. (This second peak is not due to the contribution 
of hot plasma near logr[K]?a6.2 through a Sixi line, cf. 
Fig-ffl)- The emission in Fig.|3]is calculated for looking at 
the loop apex edge-on in one AIA pixel, which corresponds 
to an integration over a height around the apex equivalent 
to 450 km, or to an arc length of «6 Mm to both sides of 
the apex. The count rate in this channel drops to very low 
values again once the condensation starts sliding down the 
loop and leaves the region of the loop contributing to the 
the AIA pixel looking at the loop apex from edge-on. As 
can be seen from Fig. [5] (bottom row), the condensation is 
at an arc length of about 54 Mm, i.e., some 6 Mm away from 
the apex at time t«143 min . This is also the time when the 
emission in the 304 A channel drops to almost zero. 

This 304 A light curve has to be taken cum grano salis, 
because the formation of the He lines in the extreme ultra- 



violet is still not fully understood, which then also applies 
to the synthesized count rates in the 304 A channel. 

4.2. Shut-off heating 

The situation is quite different for the synthesized count 
rates when the heating is shut off. The cooling sets in, and 
along with it the draining starts immediately. Thus the den- 
sity will be considerably lower in this case, leading to sig- 
nificantly lower count rates than in the the catastrophic 
cooling case. 

"Hot" channels, 193 A and 211 A. The peak in count rate 
for these channels slightly precedes the time when the tem- 
perature of maximum contribution is reached (see Fig.Q]). 
This is because the count rate is the temperature response 
weighted with the density squared and the loop is drain- 
ing quickly, leading to a rapidly falling density. More im- 
portantly, the rapidly falling density also causes the count 
rates to be very low when shutting off the heating — con- 
siderably lower than for catastrophic cooling. 
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After Smin, when the peak in 211 A is reached, the 
density has dropped already by a factor of 2. In contrast, 
during the first peak in 211 A for the catastrophic cooling 
(t?sl05niin, Fig. [3]), the density has increased by some 60% 
to 70%. Thus at the time of the (first) 211 A peak the den- 
sity in the catastrophic cooling case is higher by a factor 
of slightly more than 3 as compared to the case of shutting 
off the heating. This leads to the difference by a factor of 
about 10 in 211 A count rate during the first 211 A peak 
(compare lower right-hand panels of Figs. [3] and H]). The 
same arguments hold for the 193 A channel, which shows 
an even larger difference because it peaks later. 

Because the density drops so rapidly, the secondary 
maximum of the temperature response function at low tem- 
peratures does not play a role here. Thus, unlike the catas- 
trophic cooling there there is only a single peak in each of 
these "hot" channels. 

"Cool" channels, 131 A and 171 A. For the same reason as 
the "hot" channels these "cool" ones show much lower 
count rates than in the case of catastrophic cooling, only 
now the difference is even greater, up to a factor 1000, be- 
cause they peak later when even more material has drained. 
At this very low count rate level, they show a more com- 
plex light curve, because of the multi-peak structure of the 
temperature response at higher temperatures (cf. Fig.[T]). 

He II channel, 304 A. When employing the AIA temperature 
response, this channel shows only one single peak early on. 
This is because of the side maximum of the 304 A tempera- 
ture response near logr[K]w6.2. Consequently this peak ap- 
pears almost simultaneously with the 193 A channel peak. 
Later in the evolution of the cooling loop, the density is 
simply too low to produce a noticeable signal in this chan- 
nel. 



Switching on the heating again. After switching on the heat- 
ing rate again at i»i57min, the temperature almost imme- 
diately resumes its original temperature (Sect. [321) • During 
the subsequent filling of the loop with material through 
evaporation, the count rates in all channels creep up mono- 
tonically until they reached the initial values at t=0 (not 
fully shown in Fig. 2]). In particular, no other peaks are 
found in the count rates. 



5. Comparison to observations 

One conclusion from the discussion in Sect. 3] is that a sim- 
ple shut-down of the heating rate cannot produce signifi- 
cant observable signatures above the limb that would be 
detectable with AIA/SDO, simply because the expected 
count rates would be too low. Even if one increased the den- 
sity of the initial loop (to unrealistically higher values), the 
quick draining would prevent significant count rates. This 
rules out that simply shutting off the heating rate could re- 
produce the cooling structures seen in condensations above 
the limb. 

Catastrophic cooling provides a mechanism that pro- 
duces high count rates in the AIA bands and gives complex 
(not single-peaked) light curves. In the following we com- 
pare the synthesized count rates of the catastrophic cooling 
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Fig. 5. Comparison of synthesized (top) and observed (bot- 
tom) emission of catastrophic cooling in SDO/AIA bands. 
The count rates for the AIA bands are scaled using the 
factors listed in the top panel in order to fit better in a 
single plot. Both synthesized and observed rates are scaled 
by the same factors. The synthesized count rates are for 
an edge-on observation of the loop near the apex averaged 
over 1.3 Mm («3 AIA pixels). See Sect.O 



model to a recent AIA observation by iKamio et all (1201111 
of a condensation above the limb. 

The bottom panel of Fig. [5] shows the light curves in a 
single AIA pixel of the same structure a bove the limb as al- 
ready discus sed by Kamiq et al. (|2011f ) and shown in their 



Fig. 5. While iKamio et aD (2011) show the normalized pro 



file, we now give the actual count rates (with some scaling 
so that all curves fit into the same panel). To subtract the 
background contribution, we corrected for the count rates 
found after the condensation event (in our plot the average 
after time i«320min). The light curves clearly reveal single 
peaks in the "cool" 131 A and 171 A channels. The "hot" 
channels show multi-peaked structures, with 193 A having a 
broad maximum with maybe two or three peaks, and 211 A 
a clear double peak and some indication of a third peak 
in-between, co-temporal with the 131 A and 171 A peaks. 

In the top panel of Fig.[S] we show the count rate syn- 
thesized from the model looking horizontally at the apex 
edge on. This is similar to Fig. [31 but now we have averaged 
over a region corresponding to three AIA pixels in height 
(equivalent to 1.3 Mm). This is done to get a better match 
to the actual AIA resolution, which is (predicted to be) of 
the order of 1.6" corresponding to roughly three pixels in 
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the channels used in our study (Table 8 of H oerner et al.l . 
l201lh . As for the observations in the bottom panel, we em- 
ploy a scaling so that all curves fit into the same panel. 
The scaling for the data synthesized from the model and 
the actual observation is identical. 

Because of the averaging, the second peak in the "hot" 
channels near i=120min is not visible here. This very nar- 
row peak is strongly located in space at the place where the 
condensation occurs and is thus outshone by the emission 
along the top part of the loop we average over when looking 
at the top 1.3 Mm at the loop apex (corresponding to an 
arc length of more than 10 Mm) . The structure that remains 
in the light curves synthesized from the model is a broad 
double peak in the "hot" channels 193 A and 211 A and a 
narrower single peak in the middle in the "cool" channels 
131 A and 171 A. 

Comparing the data synthesized from the model and the 
actual observations in Fig.[5j it is obvious that some general 
features match. The "hot" channels show a broad double- 
peak structure, the single peaks in the "cool" channels are 
comparably narrow and appear roughly in the middle of the 
"hot" channel light curves, and finally the order of magni- 
tude of the predicted count rate matches the observations 
(within a factor of about 5) and the ratios of the different 
channels roughly match (within a factor of 2). While we do 
not discuss the Hen 304 A channel further because of the 
problems with our knowledge of the He line formation, i t 
is noticeable that the observation by iKamio et all (1201111 . 
their Fig. 5, shows a narrow peak in 304 A (co-temporal 
with the 171 A peak) followed by a broad second peak. This 
is similar to our synthesized light curve for 304 A in Fig.[3l 

All this is achieved without fine-tuning the model; i.e., 
we did not run a large number of models and just picked 
the best possible match. As outlined in Sect. 12. II the only 
real free parameters are the loop length L, the heat input 
Fb_, and the scale length of the heating rate A. However, 
for the comparison these are not free to choose: L is set 
by the height above the limb where the condensation is 
observed, and Fh is set by the requirement to reach a apex 
temperature above the formation temperature of the 211 A 
channel (and we chose logr[K]w6.5). This then also sets 
the pressure and the density in the loop. For A we had to 
choose a value so that the catastrophic cooling occurs. In 
the end, there is not that much room to play with the model 
parameters. 

Of course, there is also a significant difference between 
a model and observations. In the model the enhancement 
in the synthesized AIA light curves is restricted to some 
20 min in the "cool" channels and some 60 min in the "hot" 
channels. The corresponding time scales in the observations 
are about a factor of five longer. Because the time span of 
the co ndensation depends on th e the heating scale length 
A (e.g. iMliller et all l2003l . |2004[) . an alteration of this pa- 
rameter could lead to a better match. Also including of 
non-equilibrium ionization could work into this direction, 
as it will lead to lo nger time scales for the condensations 
l|Miiller et al.l . l2004h . 

6. Conclusions 

We ran two types of models for cooling in coronal loops 
and compared them to observations. In the heating off/on 
model we simply shut down the heating, which leads to 



cooling and draining of the loop, and then resumed the 
original heating. In the catastrophic cooling model we kept 
the total heat input constant but concentrate the heating 
more towards the footpoints. This leads to a thermal non- 
equilibrium near the apex, and a condensation forms in 
a runaway process. While these two processes have been 
studied in the past, we put them here in the new context of 
observations at the solar limb in extreme ultraviolet pass- 
bands by AIA/SDO. 

We found that the heating off/on model will not gener- 
ate significant count rates to understand observed conden- 
sation events above the limb. In contrast, the catastrophic 
cooling model seems to provide a natural explanation for 
the observations of condensations, with a good match to 
the observed light curves. Not only are the observed count 
rates roughly matched, but also the single peaked nature of 
the lightcurves in "cool" channels (131 A and 171 A), and 
the multiple peaks in the "hot" channels (193 A and 211 A) 
are reproduced. 

Good arguments have been presented that the pro- 
cess of thermal non-equilibrium cannot be used as a gen- 
eral process to unders tand all properties of coronal loops 
(jKlimchuk et al.l . l2010D . Nonetheless, this process is impor- 
tant in a subvolume of the corona, e.g. where the coronal 
condensations form, and it is a valuable tool for investi- 
gating the dynamic response of the corona to the spatial 
distribution of the heat input into the corona. 

The important conclusion from our study is that, even 
in condensation events in the corona, a constant supply of 
energy is needed to keep the coronal pressure. The heating 
off /on model fails, because the support of the corona is lost 
when shutting down the heating and the material drains 
very quickly. To see condensations in the corona, one has to 
keep the energy input at a sufficient magnitude to support 
the hot plasma further. In the catastrophic cooling model, 
this is achieved by a higher concentration of the heating 
towards the loop footpoints, which reduces the heat input 
in the high corona but keeps the total amount of energy 
supply to the corona at a (more or less) constant level. 
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